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I. INTRODUCTION 

Integrable quantum field theories in 1 + 1 space-time dimensions are fundamentally char- 
acterized by their factorizable S-matrix |l[|. The S-matrix along with certain bootstrap 
axioms in principle characterize the matrix elements of fields in the space of asymptotic 
multi-particle states (Form Factors) 0Q]. Finally, the correlation functions at zero temper- 
ature can be characterized by infinite integral representations using the form factors. The 
S-matrix also governs the thermodynamics of the system through the Thermodynamic Bethe 
Ansatz (TBA) 0. 

For many applications, especially in Condensed Matter Physics (see, for instance 0), 
the finite temperature correlation functions are of importance. In the Matsubara imaginary 
time formalism |0], finite temperature correlation functions can be viewed as correlation 
functions on an infinite cylindrical geometry, where the spacial coordinate — oo < a; < oo 
runs along the length of the cylinder and the euclidean time lives on a circle of radius 
R = 1/T, where T is the temperature. In a picture where the hamiltonian evolution is along 
the circumference of the cylinder, let us call it the /^-channel, the Hilbert space lives on the 
infinite line of the coordinate x, thus the states of the Hilbert space one has to sum over 
in performing thermodynamic averages have the usual multi-particle description in infinite 
volume. Deceptively simple as it may appear, this is an important observation, since it 
implies that in this picture, the matrix elements of operators do not depend on R, and can 
be expressed then by the usual, temperature-independent form factors. Formally, one has 

{0)r = ^Tr (e-^^O) (R - channel) (1.1) 

where Z is the partition function. In this paper we use the i?-channel to formulate the finite 
temperature correlation functions, the main goal being to use the Thermodynamic Bethe 
Ansatz to make sense of ( |1 . 1| ) . We show that the correlation functions can be characterized 
using only the usual form factors and some thermodynamic data, the so-called filling frac- 
tions, which are available from the TBA. We present explicit integral representations for the 
one and two point correlation functions. 



Our formulation should be contrasted to the orthogonal picture, where hamiltonian evo- 
lution is viewed as along the length of the cylinder (x-direction). In this "L-channel", the 
Hilbert space is on the finite volume R, and the computation of correlation functions does 
not involve the sum over thermodynamic states but rather involves the notion of the ground 
state \0r) on the circle of radius R 

{0)r = {0r\O\0r) (L - channel) (1.2) 

Presently, the precise structure of this ground state has not been identified and moreover, in 
this L-channel, matrix elements of operators may depend on i? in a complicated unknown 
way. There has been however some progress in this direction by Smirnov. In fact, in a series 
of papers 0, he studied properties of such matrix elements in a semi-classical limit. 

II. ONE-POINT FUNCTIONS 

Since the main features of the finite temperature calculation of correlation functions are 
already present in the simplest case of the one-point function of a quantum field operator 
O, let us start the discussion from the analysis of this case. Throughout this paper we let 
R = 1/T, where T is the temperature. To simplify the following discussion, we assume that 
in the field theory in question there is one kind of particle A of mass m. Generalization to 
other cases is quite easy. Multi-particle states are denoted as \9i,- ■■,9n), where 9 is the 
rapidity parameterizing the energy and momentum 

e{9)=mch9, k{9)=msh9. (2.1) 

Using the resolution of the identity, 

the one-point function of a local field 0{x,t) at finite temperature has the formal represen- 
tation 



((^(^,t))^ = ^E;^/^---|^ n^^^^^'^^ (^n,---,^i|0(O,O)|^,,---,^.) , (2.3) 



1 ~ 1 r dOi d9„ 

where Z = Trfe^^^j is the partition function. From translation invariance, the above 
quantity does not depend on x, t and is only a function of the scaling variable mR. 

Our aim is to show that the previous expression can actually be organized in a way which 
better reveals its physical content and which moreover provides an efficient method for its 
actual computation. In addition, we will also show that the final form of the one-point 
functions may be regarded as a generalization of the formula obtained in the case of a free 
theory, once an appropriate dictionary between the interacting and free theories has been 
established. For this reason, let us initially consider the case of a free fermionic modelQ, i.e. 
an integrable model with S-matrix S = —1. This case was originally studied in |^. For the 
following considerations, we need in particular to recall the main findings reported in the 
appendix A of 0. 

A. Main results in the free case 

In carrying out the sum over states it is important to include (5-function contributions 
to the form factors as follows. Let A = {On ■ • • ^i} and B = {6[ ■■■6m} denote some ordered 
sets of rapidities. Then the form factor is expressed as a sum over all ways of breaking up 
A and B into two sets: 

{A\Oix,t)\B) = Y. 5A^,5e,e,(A|^2)(A|0(x,t)|i3l)conn , (2.4) 

A=AiUA2;l3=BiUB2 

where S'_4^^^ are S-matrix factors required to bring |^) into the order 1^25-^1)5 i-e. (^| = 
SyX^yXj^{A2,Ai\, and similarly for \B). The inner products (^21-62) are most easily evaluated 
by introducing free particle creation-annihilation operators |^i ■ ■ ■ On) = ^^(6'i) ■ ■ ■ v4"''(6'„)|0), 



^In this paper we consider only "fermionic" theories with S{9 = 0) = — 1. The formal extension 
to bosonic theories with S{9 = 0) = 1 is a simple exercise. 



and using {A{9) , A^ {9')} = 2ti5{6 — 6'); they are sums of products of functions. The 
"connected" piece of the form factor {Ai\0{x,t)\Bi)conn is defined to be the form factor 
with no overlap, i.e. (^i|i3i) = 0. The crossing relation is thus valid and one can define 

(2.5) 

where FP in front of the expression means taking its finite part, i.e. terms proportional 
to {l/i]iy, where p is some positive power, and also terms proportional to rji/rij, i ^ j are 
discarded in taking the limit. With this prescription the resulting expression is independent 
of the way in which the above limits are taken, and it is therefore the only quantity which 
has an unambiguous physical meaning. This way of taking the limit has been already used 
in the literature (see |]10|). 

There are two effects of the 5- function terms coming from (^21-^2) when A ^ Ai. The 
first effect is due to terms involving [6{6 — d')]"^; these give rise to an overall factor which 
sums up to the partition function Z, and therefore cancels out from the final result. The 
other effect of the 5- functions is to give rise to a "filling- fraction" f{6) for integrations over 
rapidity. The final expression of the one-point functions is then given by 

(0(x,t))« = ^i-/^---^ n/(^Oe-^''M(^n---^i|0(O)|^i---^.U^ (2.6) 

where here e = e{6)R, the connected form factor is defined in (|2.5| ), and 

/(") = rrhm ■ (2-^) 

In the following it will be convenient to introduce the functions 



so that f = f^ and fe ^ = /_ . As we will show below, the expression ( p.6| ) for the one-point 
function is the one which generalizes to the interacting case. 



B. Interacting case. Quasi particle excitations 

Let us now turn to the interacting case. Again, for simplicity, we assume the theory has 
a spectrum given by a single particle A of mass m with an S-matrix S{6). We define 

a{e) = -zhgs{e), ^{6) = -z-^\ogS{e) . (2.9) 

The partition function at a finite temperature T and on a volume L (for L —>■ oo) is deter- 
mined by means of the the Thermodynamic Bethe Ansatz equations as follows [^Q . In a 
box of large volume L, < x < L, the quantization condition of the momenta is given by 
^ik(ei)L Y[-^^ S{9i — 9j) = 1, which can be equivalently expressed as 

ml sh9i + Y, ^(^i - %) = Svrrii , (2.10) 

where n, are integers. Introducing a density of occupied states per unit volume pi (6) as well 
as a density of levels p{0), in the thermodynamic limit eq. ( p.lO|) becomes 



2np = e + 27iip* pi , (2.11) 

where (/ * g){e) = /^ de'f{e - e')g{e')/2n. Defining the pseudo-energy e{e) as 

- = TXT ' (2-12) 

p 1 + e^ 

the minimization of the free-energy with respect to the densities of states leads to the integral 

equation 

e = eR-ip*log{l + e-') , (2.13) 

and the partition function is then given by 

de 



Z{L, R) = exp 



ml I —cW log (l + e-"(^)) 



(2.14) 



The interesting point is that the above partition function can be interpreted as one of a 
free gas of fermionic particles but with energy given by e{6)/R. Namely, there is a one-to- 
one correspondence between the above expression ( |2.14| ) and a partition function computed 
according to the following sum 
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oo -| „ jn JQ n 

^(L^R^ = T.^J^---^{0n---9i\9i---en) Ue-''^' , (2.15) 



1 rdOi de„ 
n=o^-J 27r 27r .^, 

where the scalar products of the states are computed by applying the standard free fermionic 
rules. To see this, let us define 

— ch^ log (l + e-^(^)) , (2.16) 

which admits the series expansion 

F{R) = j:^-^In{R) , (2.17) 

where 

In{R) = J^chee~"<'^ . (2.18) 



n=0 ^ 



de 



The partition function obtained by the TBA equation has the following expansion in powers 
of {niL) 

Z{L, R) = l + {mL)F{R) + ^^^{F{R)f + ■ ■ ■ ^-^[F{R)r + ■■■ (2.19) 



Let us now compute the partition function by using the other expression ( |2.15| ), with a 



regularization of the squares of the 5-functions appearing in {On ■ ■ ■6i\9i- ■ ■ On) as for a free 
fermionic theory on a sufficiently large volume L: 

[5{e - e')f = — ch(^) 6{e - e') . (2.20) 

Hence 

Z{L,R) = 1 + Z, + Z2 + ---Zn + --- (2.21) 

where for the first terms we have 

de ^,„ ^r.^ r de 
< £ 

de 



z, = J—<e\e> e-^(^) = J —de'5{e - e') <e'\e> e-^(^) = (2.22) 

/de 
— ch(^)e-"(') = {ml) h ; 
27r 



Z J zn zn z J zn zn <- 



-6{0i - 62)6(62 - 6^))] e"^(^^)-^(^^) = UmLYlf 



1 



(mL)/2 , 



and for the next few 



(mL)3 (rriLy (mL) 
"-'1 7: — -'i-'2 H ;^ — -fa 



3! 



(2.23) 



(2.24) 



''"^) -/J - (™l)»/?/. + (^"'^ 



I'^'^^l 



[mL) 



(2.25) 



4! "^ """"' '^"' ■ 2 

It is not difficult to carry on to liigher order and to show that the series ( |2.15| ) precisely 
coincides with the one of eq. ( p.l9| ), the only difference being in the arrangement of the single 
terms. For instance, collecting the terms proportional to [mL) in all Z„, one simply obtains 
for their sum F{R), whereas the sum of the terms proportional to {mL)"' appearing in all 
Zn gives rise to the higher power (F(i?))"'. 

The above remarkable fact can be interpreted as meaning that all physical properties 
of the system can be extracted by the quasi-particle excitations above the TBA thermal 
ground state. These excitations have dressed energy e = e{6)/ R and dressed momenta k{6): 



e{6)/R, k{6) = k{6) + 27i{(T*pi)(6) 



(2.26) 



In this context, the rapidity 6 plays then the role of a variable which simply parameterizes 
the dispersion relation of the quasi-particle excitations. This result was already understood 
in a non-relativistic situation by Yang and Yang [Q and a generalization of their result in 
our relativistic context is given for convenience in appendix A. 

In the light of these considerations, let us return now to the sum ( p.3[ ). It is now 
clear that if we replace e{6) by e{6) and proceed to evaluate the sum as in the free theory 
using ( p.4{ ), the partition function will properly be factored out and be canceled by the Z 
in the denominator, and furthermore the filling fractions will again be generated with e 
now given by the TBA pseudo-energy. Thus, one obtains precisely the same formulae as 



]6|, |2.7| ). We emphasize that the connected form factor appearing in ( p.6|) is the usual 



(temperature-independent) form factor. Then the interpretation of eq. (|2.6| ) is quite clear, 
i.e. each multiparticle state contributes to the thermal sum with its T = form factor but 
each is weighted with the thermal probability of the occupation of that state, as determined 
self-consistently by the TBA equations. 

A non-trivial check of the above claim is readily available for the trace of the energy 
momentum tensor O = T^. The connected form-factors of this operator can be easily 
extracted and they depend only on the model through the function yj. For the first ones we 
have 

(^|T;|^)conn = 27rm2 ; (2.27) 

(^2,^l|T;|01,02)conn = 47rmV(^l " ^2)ch(^i - 62) , 

and inductive application of the form factor residue equations leads to 

{On ■■■9i |T;|^i ■ ■ ■ 0„),o„n = 27rm2 if {912^623) , ■ ■ ■ </^(^„_i,„) ch(0i„) (2.28) 

+ permutations (2.29) 

where 6ij = 6i — 6j. For the aim of computing the thermal trace (i.e. considered as inserted 
into integral expressions over 6i), the last expression can be taken as 

(^„---^i|T;|0i---^„)eonn = 27rm2n!y;(^i2)v'(M , ■ ■ -(/^(^n-i.n) ch(^i„) . (2.30) 

Finally, the vacuum expectation value of this operator is given by 

2 

TTTT? 

<o|T;|o>.(r;)„ = ^^— ^. (2,31) 

where Oj are the resonance angles entering the two-body scattering amplitude of the particle 
(in the more general case, one should consider the scattering amplitude of the particle with 
the lightest mass). Therefore we have 



\n=l '' 



n JO 



..=1 2vr 



<^(^12)---<^(^n-l,n)ch(^ln)l • (2.32) 

The n-th term of this series can be represented by the graph 




where the dots are associated to the weights f~.{Oi), the thick hnes between the dots to the 
functions (f{6i — ^j+i) and the dotted hne which hnks the ending dots to ch(^i„). 
On the other hand, (T/^)i? can be computed directly from the TBA Q as 

TO«-TOo = f^[i?i?(i?)] , (2.33) 

where E{R) = — log Z/L. This can be expressed as 

(T;)^ - (T;)o = mj de-^^^ (^dR6 did - ^dee sh^) , (2.34) 

where the functions due and dee satisfy hnear integral equations which can be easily solved. 
In fact, define the kernel 

K{e,e') = y,{d,e')^^^^ . (2.35) 

Then, by differentiating the integral equation ( p.l3| ) one obtains 

(1 - K)dRe = e, (1 - K)^de6 = k , (2.36) 

where {KdRe){6) = J d9'K{9,9')dRe{6')/2TT. Introducing the resolvent L satisfying 

(l + L)(l-ir) = l , (2.37) 

we have 

Ore = (1 + L)e , ^dge = (1 + L)k . (2.38) 

Using 1 + L = J2^=o K"', one can easily express the above in an integral series: 

dRe = e + <f*(^^^e]+... (2.39) 



and similarly for dge. Substituting these series into ( |2.34] ) one then finds a perfect agreement 
with expression ( p.32| ). 
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A model with one particle whose computation of one-point functions can be explicitly 
performed is provided by the Sinh-Gordon model. This is briefly discussed in the next 
section. 

C. One point functions in the Sinh— Gordon model 



The Sinh-Gordon theory is defined by the action 
S = I £x 



1 Tfl 

{d,^f -"^ cosh g^{x) 
2 q-^ 



(2.40) 



The model is invariant under a Z2 symmetry $ -^ — $. Its two-particle S'-matrix is given 
by 

tanh \{0 — ma) 



S{e,a) 



tanh \{0 + ina) 



(2.41) 



-^. For real values of 



where a is the so-called renormalized coupling constant a{g) = 

OTT 

g the S'-matrix has no poles in the physical sheet and hence there are no bound states. The 
function (p{0) of this model is explicitly given by 

2 sin a cosh 6 



v{e) 



cosh 6 — cos^ vra 



(2.42) 



We are interested, in particular, in computing the one-point functions of the exponential 
operators of this theory, \&fc(a;) = e^^^ . To this aim, let us recall the general structure of the 



form-factors of the Sinh-Gordon model, as determined in [|11|,0 , specializing at the end to 
the operators in question. 

For scalar operators, the form-factors of the theory can be parameterized as follows 



Fn{Ol, . . . ,0n) = Hn Qn{Xl, ■ ■ ■ , Xn) _Q ', — T , 

where Xi = e^' and 9ij = 9i — 9j. -Fmin(6') is an analytic function explicitly given by 

2 



(2.43) 



^min(6',a) 



n 

fc=0 



r(fc + | + f)r(fc + i + f + f r A: + i-f + f 



r(A: + i + || r u + i-f + ii r fc + i + f + ,^ 



a I iO 



(2.44) 
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where 6 = in — 6. It satisfies the functional equations 



Fmin(^)=Fmin(-^)5(0,a) , 



i^min(^vr -6)= F^S-K + 9) , (2.45) 

^min(^7r + 6*, a)Fmin(6', a) 



sinh 6 + sinh(z7ra) 

where the last equation is particularly important for the computation at a finite temperature. 
Hn are normalization constants, which can be conveniently chosen as H2n+i = -f^i/i^", H; 



1 



,2n-2 •^.u ,, ::^ ( 4sin(7r«) \2 



2n — 

H2U^^^^, with u = ( J^'^'"}!^"\ ) ^ . In absence of additional requirements on the form-factors, 
Hi, H2 are two independent parameters. Finally, the functions Qn{xi, ■ ■ ■ , x„) are symmetric 
polynomials in the variables Xi, which are fixed by the residue equations satisfied by the form 
factors and whose total and partial degrees depend on the specific operator. In general, they 
can be expressed in terms of the elementary symmetric polynomials cr^" (xi, . . . , x„), defined 
by the generating function 



l[{x + x,) = Y.x''-'cr'j^\xi,X2,...,Xn). (2.46) 

i=l fc=0 

Conventionally the a^ with k > n and with n < are zero. It is also convenient to 

sin{n7ra) 



introduce the symbol [n] defined by [n] = HHIiZlfii.. As shown in [|I^, an interesting class 
of solutions of the recursive equations satisfied by the Form Factors (FF) is given by the 
following polynomials Qn, 

Q^{k) = \\M,,{k)\\ , (2.47) 

where Mij{k) is an {n — 1) x [n — 1) matrix with entries^ 

Mij(k)=(T2i^,[z-j + k] , (2.48) 

and the vertical lines denote the determinant of that matrix. These polynomials depend 
on an arbitrary number k and satisfy Qn{k) = (— 1)"^^(5„(— A;). Form-factors of several 



^For simplicity the dependence of Qn{k) on the variables Xi has been suppressed. 
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operators can be expressed in terms of the polynomials Qn{k). For instance, the whole set 
of FF of the elementary field $(a;) are given by (5„(0), with the choice Hi = 1/v^, H2 = 
whereas those of TNx), are given by the even polynomials Q2n{^), with the choice Hi = 0, 
H2 = 27rm2. 

The form-factors of the exponential operators "^kix) = e'^^*, are just given by the Qn{k), 
with the normalization constants fixed to be H^ = fi[k], iJ^ = ^"^[k]. Since the one-point 
function of operators which are odd under the Z2 symmetry $ -^ — $ is identically zero, 
the one-point correlators at finite temperature involving the exponential fields coincide with 
those of the even combination \E'fc + ^_fc. Using the prescription (|2.5| ) one finds 
{cosh kg^)ji ^ .r,-,o . f d9 



cosh kg^O) 
dOi de2 



l + 4[A;]2sin7ra f ^f^O) + 4[k]^ sin 



na X 



27r 27r 
dOi de2 dOs 



f-iOi)f-i02Mei2){[kYcoshei2- 



[k-i][k + i]' 

cosh 6^12 



+ 



(2.49) 



+ f l^l^l^f-i0i)f-{d2)f-{O3) (4sin7ra[fc]V(^i2)¥^(^23)cosh^i3+ 

J zvr ZTT zn ^ 

+ W(^i2,^i3,^23,A;)) + --- 

where the function W(6'i2, 9i3, 6*23) appearing in the 3-particle contribution can be expressed 

as 



2 <^ (^12) <^ (^13)75(^23) X 
1 
cosh 6^12 cosh 6*13 cosh ^23 

COsh(gl3 - ^32) + COsh(gi2 - ^23) + COsh(g21 " ^13, 

cosh 6^12 cosh 9i3 cosh 6'23 
with the constants A{k), B{k) and C{k) given by 

nk? A, ,,2r. , .12 r. .ir.i2r. , .A , i^f 



W(^12,^13,^23,A;) 

A{k) + B{k) 
+ C{k) 



(2.50) 



A{k) 



B{k) 



3roi2^ 



[k] 
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[\k - lf[k + If -[k- l][kf[k + 1]) + i^ ([A; ~2][k + If + [k- lf[k + 2] + [A;]=^[2] 
{^\k -2][k- l\[k + 1]=^ + [k- 2][kf[k + 2] + [A; - lf[k +l][k + 2]- [kf{l - [2f)) 



[k] 



16 



[[k - 2\[kf[k + If + [k- lf[kf[k + 2] + 2[k - lf[k][k + If + [kf+ 
[k - 2\[k -l][k + If -[k- 2\[kf[k + 2] - [A; - lf[k + 1][A; + 2] + 
[k - l][kf[k + 1] - 3[A; - 2][A; - 1][A;][A; + 1][A; + 2] ' 
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c{k) = -\[k-i][knk+i] 

The expressions for higher terms are rather comphcated and are not given here. However, 
one can easily check that for /c = ±1, the above expression coincides with the corresponding 
series of the trace of the stress-energy tensorQ. The plots of some of these quantities for 
some values of the coupling constant are shown in Figure 1. Notice that if the operator \&fc 
is "on resonance" with the value of the coupling constant, namely if ka = n, where n is an 
integer, its one-point function coincides with its vacuum matrix element only. Moreover, 
varying a, there may be a swapping of the profiles associated to the one-point functions of 
different \&fc- 

D. Low and high— temperature limits of the one point functions 



Let us discuss in more detail the behavior of the one-point correlators ( p.6| ) as functions 
of the temperature. First notice that if the anomalous dimension of the operator O is given 
by 2Ao, its vacuum expectation value may be expressed as (0|(9|0) = {0)q = Yrn?^° , 
with Y a pure number. Hence, the series ( p.6| ) relative to its one-point function at finite 
temperature can be cast in the scaling form 

7^ = H{mR) , (2.51) 

where H{mR) is a function of the scaling variable rriR. 

In the low-temperature limit R -^ oo, the pseudo-energy goes as e{6) ~ mRcosh.6 and 
therefore the A^-th term of the series entering H{mR) vanishes asymptotically as e"^"^-^, so 
that the leading correction is given by 

H{mR)^l + ^^^^^^Ko{mR) + --- (2.52) 

where Kq{x) is the usual Bessel function, and \A) = \9). 



^For the check of the 3~particle contribution one needs the trigonometric identity [3] = —1 + [2]^ 
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On the other hand, in the high-temperature hmit R ^ 0, the one-point correlators may 
become scahng invariant functions^ i.e. 

{0)r - -^ , (2.53) 

with a power-law exponent t] ruled by the underlying Conformal Field Theory. This is the 
case, for instance, of massive integrable models obtained as a deformation of a Conformal 
Field Theory by a strong relevant field. Under this hypothesis, the limit i? — > can be 
controlled by means of a Conformal Perturbation Theory on the cylinder. Let 

A = AcFT + ^ j ^{x)d^x , (2.54) 

be the action of the off-critical model, where $(a;) is the relevant operator which gives rise to 
the massive integrable theory, with A ~ -m^^^^*. The general structure of the perturbation 
series for the one-point functions is then given by 

n 

{0)r =J2dn , (2.55) 

n=0 



where 



dn=^-^ I (e|0(O)$(Xi)---$(X„)|e)connCi'Xi---t/2X„ , (2.56) 

n\ Jcvi 



and Xi are points on the cylinder and the connected correlation functions are calculated in 
the unperturbed CFT. The field ^ entering eq. ( |2.56| ) is relative to the conformal operator of 
lowest anomalous dimension in the theory which plays the role of the vacuum state on the 
cylinder (for a unitary model, this field is the identity operator). Notice that the first term 
of the series is given by 

{^\O\0 = (^) Qos , (2.57) 



^This is not always the case, because there are theories, hke the thermal Ising model, with 
logarithmic behavior or models, like the Sinh-Gordon one, which present logarithmic corrections 
which may spoil the pure power law behavior given in the text. 
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and this will be different from zero only if the conformal structure constant C^o^ does not 
vanish. In this case, the exponent rj in ( |2.53| ) will coincide with the anomalous dimension of 
the field O itself, rj = 2Ao, otherwise the exponent t] will be given hj rj = 2Ac) — g(2 — 2A<j,), 
where the integer q is the first non-vanishing term in the series ( |2.55| ). 

The above considerations open up the possibility to extract conformal data out of the 
high-temperature limit of the one-point function. This is easily checked by the analysis of 
some simple models which have only one massive particle in the spectrum. 

1. Yang-Lee model 



m 



The S'-matrix of the perturbed Yang-Lee model was determined in |I3| to be 

tanhi(6' + 2^) 
S{d) = ..„JL JJ, ■ (2.58) 

Therefore for '~p{9) we have 

^w = - _^:;7""i • (2.59) 



v3 cosh 6 



cosh^^-i 

The only off-critical primary operator of the theory coincides with the trace of the stress- 
energy tensor. In the conformal limit it goes to the field with the lowest anomalous dimension 
2A = —2/5 which simultaneously plays the role of the vacuum of the cylinder. Hence, in 
this case we expect to find for the exponent r] = —2/5. In order to extract this parameter, 
by taking the logarithm of both terms in ( p.53| ) 



— 7ylogi?= \iia{0)ji 



(2.60) 



7] is easily identified by isolating the terms proportional to — log R in the series of the right 
hand side. Observe that, factorizing the term (T^)o in eq. ( p.32| ), the logarithm of the series 
( p.32| ) can be expressed in terms of the usual cluster expansion 



97r Tn 
log(T;)« - log(T;)o = j;^ 



d0 

2^ 



fid)e 



-e{9) 



dOi 02 



X 



ifie^ - 62) cosh(^i - 62 



27r 27r 
1 /27rm2 



f{di)f{02)e 



-e(9i)-e(e2) 



X 



2 UTl' 



A^JO, 



+ 



(2.61) 
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where, for this model 

27rm^ 



.^M^ -2v^ • (2-62) 

Following the TBA analysis [^, it is known that in the limit R -^ 0, the pseudo-energy e{6) 
tends to a plateau of value cq = log ^^^ in a region \9\ < 21og(2/i?). Therefore the first 
estimate of r] is obtained by the first term of the right hand side of ( |2.61|) 

(i/Tjo^ 27r (i/r)o J-iog(^) 

27rm2 1 1 , / 2 
log 



(T/r)o 1 + e^o vr ^ Vi? 



I.e. 



77 ~ ^^ = -0.421178... (2.64) 

which is within 5% to the actual value t] = —0.4. The next correction is obtained by isolating 
the term proportional to — log-R in the next term of the cluster expansion ( |2.61| ), which is 
equivalent to computing the quantity 

3 / 1 \^ /■+°° 1 

— ^ / de . = 0.0268126... (2.65) 

vr2 \3 + V5j J-oo cosh^e-i ^ ^ 

so that the estimate of rj improves to the values 

ri ~ -0.394365... (2.66) 

which differs from the actual value by 1%. 

The main result of these computations can be summarized by saying that the exponent 
rj will be expressed by a series (of alternating signs, for this model) whose terms involve, in 
addition to the ratio ( |2.62| ), inverse powers of vr and the plateau value eo. 

2. The non-unitary model M.3,5 

Another integrable model with only one massive excitation to which we can easily apply 
the above consideration is the off-critical "thermal" deformation of the minimal model 

17 



^^3,5) which has been previously studied in [p!5| -[T7|. The S'-matrix of the model is given by 



S{9) = — itanhi (6' — if ) and therefore v'(^) = 1/cosh^. In the conformal limit, the trace 
of the stress-energy tensor tends to the conformal field O of anomalous dimension 2A = 2/5 
and this field has a non-vanishing structure constant with the field ^ of anomalous dimension 
2A^ = —1/10 which plays the role of vacuum state of the cylinder. Hence we expect that the 
exponent t] of this theory should be equal to 2/5. By using the same arguments as above, 
with the additional data 

"v^+l" 



27rm^ , 

TTTvir = 2 , eo = - log 

for the first estimate of the exponent rj we have 



(2.67) 



n ~ - ^ V- = 0.393453... (2.68) 

vr3 + y5 

which is within 1% of the actual value t] = 0.4. In this model, the fact that the first term of 
the cluster series gives a better approximation to rj seems to be related to the vanishing of 
the next leading correction in the cluster expansion. 

III. TWO-POINT FUNCTIONS 

Having understood the main features of the one-point functions, the same reasoning 
applies to the higher point functions. In this section we focus on the two-point functions. 
We begin with the double summation over states 

{o(x,t)omn = ^T.e''''''ii^p{x,t)mmo{Qm . (3.1) 

As explained in , for a free theory the terms in ( p.4|) with Ai ^ A give rise to three kinds 

of contributions: 

(i). Terms that diverge in infinite volume involving [6{6 — 9')]'^ which can be arranged into 

an overall factor of Z. 

(ii). Terms that give rise to factors of the one-point function {{0)rY. 



(iii). Terms that modify the integration over rapidity for both the ip and ip' states by the 
same fiUing fraction /. 

The remaining contributions are only expressed in terms of connected form-factors for which 
crossing symmetry is vahd. As for the one-point function, the only difference between the 
free and interacting theory is that e, k should be replaced by the TBA dressed quantities 
e, k. A further justification of this will be given below. 

Introduce an index a = 1 for the ip' particles and a = —1 for the ip particles. For the ip' 
particles the integrations J d6 are accompanied by f{6) defined in (|2.7|) , whereas for the i/j 
particles one has f(9)e~'^^^^ due to the e~^^^ in (|3.1| ). These factors can be expressed as f±{0) 
as defined in ( p.8| ). By using the crossing symmetry, the connected form factors appearing 
in ( p.l|) are ('?/''|O(0)|'?/')conn = ('^ + ^7r,'?/^'|C(0)|0), where ip + m denotes all rapidities shifted 
by m. We can describe this by defining 

(0|O|^i ■ • ■ 9^).,....^ = {0\O\e, - inai, ■■■6^,- ina^) , (3.2) 

where a = (a - l)/2 G {0, 1}. Finally, 



1 ^ fdOi dON ^ 



io(.,t)om.-m.fZj^_Zji^ 









(3.3) 



where e^, kj = e{9j), k{6j). 



A non-trivial check of the above expression is the Kubo-Martin-Schwinger formula pO 
which for a bosonic operator O reads 

(C»(x, t + R)OmR = (0(O)0(a;, t))^ . (3.4) 

Shifting t by i? in ( |3.3| ) gives the factor e"'^'^. Noting that fae~'^^ = f-a, we simply relabel 
the particles a ^ -a in the sum. Since {0{Q)0{x,t)) = (C(-x, -t)O(O)), then (^) holds 
as long as 

|(0|O(0)|^i ■ • ■ 0^).,...,!' = mom9^ ■ ■ ■ 9n)_^,..._^J . (3.5) 
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The latter identity can be proven as follows. Let | 9 — ina) denote 16*1 — inai, ■ ■ ■ ,6]^ — 
i7ia]\f) and ( 9 + i7ia\ denote {9^ + iiro-N, ■ ■ ■ ,9i + inail. Then, 

\{Q\O{0)\9i---9N).,...aJ = {0\O\t-tna)C9+t7Ta\O\0) . (3.6) 

Under the transformation a -^ —a, one has a -^ —a — 1. For local operators, the form 
factor is invariant under a shift of all rapidities by the same constant, thus 

|(0|O(0)|^i ■ • ■ ^-.,...-.^1' = {0\O\t + nxa)C9 - i7^a\0\()) (3.7) 

= {Q\0\t -ma){t + nid\0\Q) , 

where we have used crossing symmetry in the second line. The orders of rapidities in ( p.7| ) 
can be brought to the order in (gj) using (0|O| ■ ■■9i,9j, ■ ■ ■) = S{9i-9j){()\0\ ■ ■■9j,9i, ■ ■ ■). 
In doing this one always encounters the pair of factors S{9 — m)S{—9 — m). Using crossing 
S{9) = S{i7T — 9) and unitarity S{9)S{—9) = 1, the latter factor equals 1. Thus, ( |3.5| ) holds. 

IV. CONCLUSIONS 

We have proposed explicit integral representations for the one and two-point correla- 
tion functions at finite temperature involving form factors and the TBA pseudo-energy- 
momentum. Several checks were performed, in particular we showed that for the trace of 
the energy momentum tensor, our expression for the one-point function coincides precisely 
with the TBA result. We also checked that the Kubo-Martin-Schwinger relation holds for 
the two-point function. 

We close by listing some possible applications of our formalism that deserve further 
investigation. By extending our results to boundary field theories at finite temperature, 
as was done for the free case in [^, it should be possible to study the finite temperature 
properties of conductivities in quantum wires with impurities. Another application is to the 
study of finite temperature crossovers in the vicinity of a zero temperature quantum phase 
transition (for a review see 0]). Finally there is the issue of spin-diffusion in quantum spin 
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chains, which would require the study of the infinite temperature hmit of the continuum 
version of a spin-spin two-point correlation function [pj|]. 

Acknowledgments. We would like to thank F. Lesage for discussions. We are also grateful 
to P. Simon for his numerical support in producing Figure 1. 
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Appendix A. TBA Dressed Energy and Momentum 



Let {uj, 6j) and (ra^-, 6'^) be two Bethe states both satisfying ( p.lO|) , where n'^ = rij except 
for j = a. Subtracting the two equations one obtains 

mL (sh^^. - sh^j) = J2 ^(^j - ^d - ^(^j - ^D • (4-1) 

i 

{j 7^ a). Since 6'' ^ 6j we can introduce a function x(6') and write 

L (sh^;. - sh^j) ^ x(%)ch% • (4.2) 

In this thermodynamic hmit, ( [4.1| ) can be written as 

27r(l - K)ipx) = aiO - 9^) - a(0 - 9'J , (4.3) 

where p is the previously introduced density of levels and K is defined in ( ^.351 ). One also 
has 

a{e - e^) - a{e - e'J = T" d9'K{e, ^')(1 + e^^''^) . (4.4) 

J On 



Using the resolvent operator, as defined in (|2.371 ), one has 

"' dO' 



PX{0)= ["^L{0,e'){l + e<''^) . (4.5) 

Jda 277 

Now consider the difference in energy AE between the two Bethe states: 

AE = m ch^' - m ch^„ + m ( dO sM ^^^>P^^J . (4.6) 

Substituting (U), and using the property L(6', 6'')(1 + e^(^')) = (1 + e^(^))L(6'', 6*), along with 
( p.38| ) one finds 

^E = ^{e{e'^) ~ e{e^)) . (4.7) 

The dressed momentum proceeds similarly: 

AP = Y,m sh^^. - m sh9j = (4.8) 

X{0)p{0) 



/v- 
dO ch.6 — 



22 



Again substituting ([4 .51) and using (|2.38|) one finds 

AP = k{9'J - k{9^) (4.9) 

where k is defined in ( |2.26| ). 
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Figure la. One-point correlators of e^^* as functions of mR for a = 1/20. 




Figure lb. One-point correlators of e'^^* as functions of mR for a = l/v7- 
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